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A recently formulated statistical mechanics method is used to study the phase transition occurring 

r -1 " in a generalisation of the Traveling Salesman Problem (TSP) known as the centred TSP. The method 

is^ , shows that the problem has clear signs of a crossover, but is only able to access (unsealed) finite 

■ temperatures above the transition point. The solution of the problem using this method displays a 

' curious duality. 

(D 

00 ■ I. INTRODUCTION 

Ch I In a recent paper, Lipowski and Lipowska Q considered a generalisation of the geometric TSP by adding a new 

C ■ term to the cost function: 
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H = {l~a)HL+aHc (1) 

Hl = J2^L{v^+uV,) = Y,\r^+l-n\ (2) 



He = ^yc(r,+i,r,) = ^|r,+i+r,| (3) 

Q The familiar term of the traditional TSP is the total path length of a cyclic torn' defined by some permutation of 
Q the N city indices. The new term He is a measure of the distance of the tour from the centre of the geometry. Here 
we shall only consider the two-dimensional problem in a square domain, and shall take this centre, the origin of the 
I ■ city coordinates, r^, to be the centre of the square. The parameter a, < a < 1, specifies the relative strength of 
^ . each term. 

' The authors of ^ gave example motivations for the centred model, but the fact that they found a transition 
' as a is varied and were able to identify it as a complexity transition provides sufBcient interest. Specifically, they 
i used simulated annealing techniques to numerically study zero temperature optimum configurations. While these 
configurations appear similar to standard TSP tours for small a, the picture is rather different at large a. For a — 1, 
i; rather than connecting near neighbour cities (r^+i ~ rj), the tour connects cities lying on points that are almost 

■ inverted with respect to the origin (r^+i ^ — r^). This change in character can be seen directly though an observable 
-j~>, , measuring the angle between successive links of the tour: at small a the angle is not especially biased in any direction 

' (though boundary effects cause a slight preference for acute angles) , while at large a there is a distinct preference for 
Ch , links to double back on themselves. Sharp changes in quantities such as (Hl) as a is varied suggest a phase transition. 
I ■ Moreover, this transition was identified as a complexity transition since the simulated annealing schedule required 
' O ' to reach optima in the (large a) C-phase is much less stringent than that required in the (small a) L-phase, indicating 
^ . a difference in characteristic difficulty of the problem as the transition is crossed. Complexity transitions have recently 
been of strong interest to computer scientists, combinatorial optimists and physicists, and although a transition has 
already been identified |3| in other versions of the TSP, the geometric basis and simplicity of this centred model are 
very appealing. 

Many questions about the nature of the transition remain, and the aim of this paper is to investigate the model using 

■ analytic techniques recently developed to study the statistical mechanics of TSP-like models These techniques 
[ allow a full solution of the theory in the region of high city density and finite temperature. In particular, expectation 

values for the path length (Hl) and its fluctuations can be computed and these provide evidence for a smoothed 
transition at finite temperature. The technique is not able to access the optima at low temperature and we cannot 
use it to compute critical indices for example. Nonetheless, our methods produce other insights such as a duality; 
a <-> (1 — a) between thermodynamic quantities that may have a role beyond the region of applicability of the 
technique. 
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FIG. 1: The average link length (l/N) (Hl) for the centred TSP according to the formalism of ^ shown for (3 = 2.0, 5.0, 10.0, 20.0 
with increasingly sharp step. Points (omitted at /3 = 20.0) are from Monte Carlo simulations. 



II. FORMALISM 



The analytic techniques we use to solve the model were derived in using a functional formalism and have since 
been placed on a firm basis via a discrete approach Q . The continuum equations for a system of unit area and a flat 
distribution of cities allow the free energy to be written in terms of a quantity s(r) which is defined by an integral 
equation, 



(3F 
s{r) 



2N / (fr logs(r) -NlogN 
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r ^ exp {-13 [(1 - a)14(r, r') + aVc{r, r')]) 



(4) 
(5) 



Observables can now be obtained by standard thermodynamic relations. For example the expectation values for the 
total length of the path and its fluctuations are given by. 



(Hl) = N / d^rd^r 



exp {-(3 [(1 - a)VL{r, r') + aVc{r, r')]) 
(r)s(r') 



a d{HL) d{HL) 
(3 da dp 



(6) 
(7) 



These expressions are proportional to N in contrast to the expected dependence of the optimum configurations. For 
example, the optimum path length for the standard {a — 0) TSP grows as \/N, and for a > has more complex N 
dependence reported in 0. Along with the non-extensive A^logA^ entropy, this scaling makes it already clear that 
the region within which the technique is valid is restricted. Similar formulae hold for expectations of He, but in all 
cases we find that the solution of the integral equation , obeys a duality relation: 



si„Q(r) = s„(r) = Sa{-r) 



(8) 



This duality is a consequence of the symmetry that inverts with respect to the centre, r — » — r, and interchanges 
Hl ^ He- As a result, F{1 — a) = F{a) so expectations of He are simply related to those of Hl, for example, 
{Hc{a)) — {Hl{^ ~ ct))- One immediate result of the duality is that if any transition occurs, it must be at Oc — 1/2. 
We shall return to discuss whether this duality holds outside the range of validity of the equations ((SJ . 

It has not been possible to analytically solve the integral equation ij^l, but an iterative procedure on a discretised 
version converges rapidly and is numerically stable. Monte Carlo simulations independently confirm the results and 
are shown in the figures below. At the lowest temperatures it becomes hard to equilibrate the Monte Carlo and we 
omit the (3 = 20.0 data. In figure Q the expectation value of the length per link is shown as a function of a, for 
a variety of temperatures. Evidently this tends towards a step function and is a smoothed out version of the zero 
temperature result shown in reference 0. 
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FIG. 2: The average per link length fluctuations {1/N){{Hl — {Hl))"^) for the centred TSP according to the formalism of 
3 shown for 13 = 2.0, 5.0, 10.0, 20.0 with increasingly pronounced peak. Points (omitted at /3 = 20.0) are from Monte Carlo 
simulations. 



Fluctuations in the path length also grow as N and figure 10) shows how {1/N){{Hl — {Hl))"^) changes with a. 
This quantity is numerically more delicate than the simple average, but although it can be computed by solving linear 
equations for ds/d(3 and ds/da, we have simply evaluated the expression in iQ with discrete derivatives. The Monte 
Carlo also suffers from large sample to sample fluctuations at low temperature. As the temperature is reduced, a peak 
starts to appear at Oc. 

Other quantities such as correlations between angles between successive links have also been computed within this 
formalism and agree with Monte Carlo simulations. The computations are interesting in that they show how to 
generalise the methods of 4] to compute correlations along the path, however they do not throw any further light on 
this problem and its transition. 



III. LOW TEMPERATURE LIMIT 



Although no general analytic solution of the integral equations has been found, the low temperature limit can be 
understood using a saddle point at large (3. Because the formalism does not scale temperature 6] this low temperature 
limit does not correspond to the TSP optimum. 

Writing Saiic) — ta{r)e~^^''^'^\ the stationarity condition of the exponent at large (3 is: 

Wa{r) = mirir' [(1 ~ a)\r — r'\ + a\r + r'\ — ^^(r')] (9) 

A solution exists with Wct(r) linear in r, matching the tendency of the iterative numerical solutions at low temperature. 
The prefactor tair), can also be determined based on linear rather than the usual quadratic fluctuations 4]. The 
resulting solution obeys duality: 



exp (-/3ar) < a < 1/2 
^ exp (-/3(1 - a)r) 1/2 < « < 1 



The resulting expectation values have sharp changes at the transition point. For example, 

1.^ r ^exp(-/3«r) 0<a<l/2 

A^^ l2/d2r|r|^(V2 + log(l + V2))/3 l/2<a<l ^ ' 

The length fluctuations develop a delta function singularity at ckc, with leading coefficient proportional to 1//3, 

The height of the fluctuation peak as determined by the iterative solution of the integral equation at fixed ac appears 
first to grow as /? as the temperature is reduced, but eventually this 1//3 behaviour is observed. 



IV. CONCLUSIONS 



The method of Q for studying the statistical mechanics of TSP-Hke problems has been applied to the centred TSP. 
Both numerical solutions of the equations at finite temperature and the limiting behaviour at large (3 provide evidence 
supporting the transition observed though numerical study of optimum configurations in Q . The method is not able 
to reach the region where optimal configurations dominate and the character of the evidence is a smoothed signal of 
the transition indicating that the transition point lies at a temperature below that accessible by this technique. 

A noteworthy feature of the solutions is a duality, F{\ — a) = F{a), relating observables on either side of the 
transition. This duality can be understood in terms of the density of states: the number of cyclic tours with Hl 
between E and E + dE is the same as the number of tours with He between E and E + dE, where in this regime 
of high city density and high temperature, E should not be too small. This kind of relationship also exists between 
a particular pair of one-dimensional Hamiltonians Q in which the configurations are also specified by permutations 
of N indices. However, in that case, even at finite N, a map between a configuration with a certain value of one 
Hamiltonian and another configuration with the same value for the other Hamiltonian was demonstrated. A brief 
consideration of particular instances of the TV = 4 case shows that this precise relationship cannot hold for the centred 
two-dimensional TSP, so the duality must only be valid in the large N limit. Within the large N limit, even if the 
duality holds for optima and near optima outside the regime of high temperature, the lack of a clear map prevents 
one from using the a ^ 1 model on the easy side of the transition to solve the hard TSP at a = 0. 
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